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The thermodynamic parameters of the superconducting state in PtH at the pressure 76 GPa 
have been examined. The calculations have been carried out in the framework of the Eliashberg 
formalism. It has been found that the critical temperature (Tc) changes in the range from 30.6 
K to 16.8 K, depending on the assumed value of the Coulomb pseudopotential: /i* £ (0.1,0.3). 
The other thermodynamic quantities differ significantly from the predictions of the classical BCS 
theory. In particular: (i) The parameter 2A (0) /fcsTc reaches the values from 4.25 to 3.98, where 
A (0) denotes the low-temperature order parameter, (ii) The ratio of the specific heat jump (AC) 
to the specific heat in the normal state (C N ) takes the values from 2.07 to 1.97. (iii) Finally, 
T C C N (T c ) /H c (0) G (0.145, 0.156), where H c (0) is the low-temperature thermodynamic critical 
field. 
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I. INTRODUCTION 



The electron-phonon interaction can lead to the formation of the superconducting state that is characterized by a 
high value of the critical temperature (Tc). In the first case, it is relevant for the multi-band systems, in the second 
case, for the high-temperature superconductors (cuprates), and in the third case, for the materials under the influence 
of high pressure (p) . 

An example of the multi-band system, in which the phonon-mediated superconducting state is induced at a high 
value of the critical temperature, is magnesium diboride [l[. On the basis of the conducted experiments, it has been 
stated that pc] MgB2 equals 39.4 K. From the physical point of view, the anomalous thermodynamic properties of the 
considered superconducting state are related to the existence of two energy gaps; the first one value (7.0 — 7.4) meV 
corresponds to the two-dimensional a band, the second one (2.6 — 2.8) meV exists in the three-dimensional ir band 

In the case of cuprates, some authors have assumed that the electron-phonon interaction, supplemented by the 
electron-electron-phonon interaction, can be responsible for the creation of the superconducting phase Q, Q, 0]. 
Note that cuprates are the materials characterized by the highest experimentally verified values of Tc- For example, 
for the compound HgBa 2 Ca2Cu308+ v (Hgl223), the following has been obtained Tc — 135 K wherein, under the 
influence of the pressure (p ~ 31 GPa), the critical temperature in Hgl223 increases to the value of about 164 K 

As it turns out, the pressure may be the factor that causes the strong increase of the electron-phonon interaction 
[hJ. In particular, it is clearly visible in the case of lithium and calcium, where the maximum value of Tc is equal to 
14 K (p = 30.2 GPa) and to 25 K (p = 161 GPa), respectively [ll|, [3. 

Additionally, it is believed that a very high value of the critical temperature can be obtained in metallic hydrogen 
[i~3l |. For instance, for the pressure at 450 GPa, the superconducting state in the molecular phase of hydrogen is 
characterized by Tc — 240 K [l4j]- At the extremely high pressure (p ~ 2 TPa), the critical temperature may be 
even equal to 600 K [H, 0. It is worth noting that from the physical point of view, the high value of the critical 
temperature in hydrogen is associated with the small mass of the atomic nuclei forming the crystal lattice and the 
lack of the inner electron shells, which results in the high value of the electron-phonon coupling constant. 

In order to obtain the superconducting state with a high value of the critical temperature, but at the pressure that 
is significantly lower than the metallization pressure for hydrogen (p ~ 400 GPa), Ashcroft proposed the method of 
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the chemical pre-compression [17] . The suggestion made by Ashcroft caused a rapid increase in the number of the 
publications on this topic. In particular, the following chemical compounds have been taken into consideration: CH4, 
GeH 4 , SnH 4 , and S1H4 [H], p]|, [U, [2l|, [H, (H- We can notice that the obtained theoretical results have been 
so promising that this branch of physics is now one of the most intensively developed. 

In the presented paper, we have discussed the results obtained for PtH under the pressure at 76 GPa. In particular, 
we have determined the most important thermodynamic parameters of the superconducting state in the framework 
of the Eliashberg formalism [24| , [25[ . 



II. THE ELIASHBERG FORMALISM ON THE IMAGINARY AXIS 



A. The Eliashberg equations 

The Eliashberg formalism on the imaginary axis (i = y— 1) allows us to determine the values of the order parameter 
function (<p n = </>(iw„)) and the wave function renormalization factor (Z n = Z (ito n )). The symbol w n denotes the 
n-th Matsubara frequency: uj n = (7r//3) (2n — 1), where j3 = (fc^T) - and ks is the Boltzmann constant. We notice 
that the order parameter is defined by the ratio: A„ = (f> n /Z n . 

The functions 4> n and Z n can be calculated using the equations below [26|, [27j, (28j : 



A (iuj n - iw m ) - n*9 (uj c - \u) m \) 



- E 

P m=-M \/^li Z l 



(1) 



n ' m = — M V lu m^ J m ~ 

where the pairing kernel for the electron-phonon interaction is given by the formula: 

A(z)E2|"" ( i!]^a 2 f(!]). (3) 

In equation ([3]), the a 2 F (0) denotes the Eliashberg function. For the PtH compound at 76 GPa, this function has 
been calculated in the paper [25j. The value of the maximum phonon frequency (fi ma x) is equal to 117.81 meV. 

The electron correlations have been modeled by the Coulomb pseudopotential /j,*; 9 denotes the Heaviside unit 
function and lo c is the cut-off frequency; lo c = 3f2 max . This study refers to a wide range of the Coulomb pseudopotential: 

fi* e (0.1,0.3). 

The Eliashberg equations on the imaginary axis compose the infinite system of the equations. However, for the 
temperatures greater than the temperature of zero Kelvin one can reduce their number due to the fact that 4> n and 
Z n become saturated at large values of n. In particular, the following range of the temperature has been assumed: 
T G (T = 2.5K,Tc). In the present case, 2201 Matsubara frequencies (M = 1100) should be taken into consideration. 

The Eliashberg equations have been solved using the iterative method, presented and tested in the papers p9j and 

H. 



B. The order parameter and the wave function renormalization factor 

Figure [TJ (A) presents the dependence of the order parameter A m= i on the temperature and the Coulomb pseu- 
dopotential. The charts have been drawn up on the basis of the results presented in Figs. [TJ (B)-(D). 

It can be noted that the value of the function A m= i (T, strongly decreases with the increasing temperature and 
the Coulomb pseudopotential. As it turns out, the course of the order parameter can be parameterized using the 
simple formula: 



where T = 3.35, and: 



A m=1 (T, M *) = A m=1 (/i*h/l- ( — j , (4) 



A m=1 O*) = 33.57 (X) 2 - 26.56fi* + 7.80. (5) 
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FIG. 1: (A) The order parameter A m= i as a function of the temperature for the selected values of the Coulomb pseudopotential. 
The circles, squares and triangles correspond to the exact numerical solutions of the Eliashberg equations. The solid lines 
represent the results obtained by using analytical scheme (see Eq.Q). (B)-(D) The order parameter on the imaginary axis for 
the selected values of T and fi* . 



It should be noted that the critical temperature is also strongly dependent on fi*, and it belongs to the range from 
30.6 K to 16.8 K. 

From the physical point of view, the quantity of 2A m= i (T, fi*) with a good accuracy determines the value of the 
energy gap at the Fermi level. Thus, the formula (|S|) is relatively accurate to estimate the physical value of the 
Coulomb pseudopotential on the basis of the tunnel experiment [3l|. 

Figure^] (A) presents the dependence of the wave function renormalization factor (Z m —i) on the temperature and 
the selected values of the Coulomb pseudopotential. These charts are based on the data collected in Figs. [5] (B)-(D). 

In contrast to the order parameter, the function Z m —\ (T, fi*) grows together with the increase of T and fi* . However, 
the changes in the value of Z m= i (T,fi*) are not too large. 

It is worth noting that the wave function renormalization factor for the first Matsubara frequency can be parame- 
terized with the use of the following expression: 

Z m= i (T, n*) = [Z m=1 (T c ) - Z m =i (fi*)] 

+ Z m=1 (p*). (6) 

The value of Z m= i (Tc) can be calculated on the basis of the formula: Z m —\ (7c) = 1+A, where A denotes the electron- 
phonon coupling constant: A = 2 L max a 2 (f2) F (f2) /f2. For the PtH compound, we have obtained: A = 1.24. The 
function Z m=1 (fj,*) has the form: Z m=1 (fi*) = -0.764 (fi*) 2 + 0.622/x* + 2.058. 

From the physical point of view, the value of Z m= \ (T,fi*) with a good approximation determines the dependence 
of the electron effective mass (m*) on the temperature and the Coulomb pseudopotential: m* = Z m —\ (T,fi k )m e , 
where the symbol m e represents the electron band mass. Based on the results presented in figure [2] (A), it can be 
noticed that the effective mass of the electron in the PtH compound is high in the entire range of the considered 
temperatures. 




4 




• n*=0.1 

2.10 - ■ H*=0-2 

A * n*=0.3 

5 10 15 20 25 30 
T K 




100 200 300 400 

(o meV 




100 200 300 400 
to meV 



FIG. 2: (A) The wave function renormalization factor Zjyi — 1 cLS Qi function of the temperature for the selected values of the 
Coulomb pseudopotential. The circles, squares and triangles correspond to the exact numerical solutions of the Eliashberg 
equations. The solid lines represent the results obtained by using analytical scheme (see Eq.©). (B)-(D) The renormalization 
factor on the imaginary axis for the selected values of T and p* . 



C. The thermodynamic critical field and the specific heat 

The thermodynamic critical field and the specific heat can be determined by calculating the difference in the free 
energy between the superconducting and the normal state [HI ]: 



(V^+Al-K|) (7) 



where p (0) is the value of the electron density of states at the Fermi level. The symbols Z% and Z% denote the wave 
function renormalization factor for the superconducting state (S) and the normal state (N), respectively. 

In the lower panel of figure [3] (A), we have plotted the dependence of AF/p (0) on the temperature for the selected 
values of the Coulomb pseudopotential. It has been found that with the increase of the parameter p* , the absolute 
values of the difference in the free energy strongly decrease. The discussed trend can be characterized by the ratio: 
l^F (To)]^.! I [AF (T )]^ a3 = 3.55. 

The next step is to determine the thermodynamic critical field: 



= v/-87r[AF/p(0)]. (8) 

vp(°) 

The upper panel of the figure[3](A) shows the obtained results. Similarly as it has been observed for the free energy, 
the increase in the Coulomb pseudopotential causes a severe decline in the value of the thermodynamic critical field: 
[He (0)1^=0.! / [H C (0)]„. =0 . a = 1.89, where H c (0) H c (T ). 

As the next step, the specific heat of the normal state has been estimated: fcBp ( j = -j), where the Sommerfeld 

constant is given by: 7 = |7r 2 (1 + A). On the other hand, the specific heat of the superconducting state should be 



5 




FIG. 3: (A) The difference in the free energy (lower panel) and the thermodynamic critical field (upper panel) as a function 
of the temperature for the selected values of the Coulomb pseudopotential. (B) The specific heat of the superconducting and 
normal state as a function of the temperature for the selected values of the Coulomb pseudopotential. 



calculated using the formula: C s = C N + AC, where the difference between the specific heat of the superconducting 
and the normal state (AC) equals: 

AC(T) = l rf 2 [AF/p(0)] 

k BP (0) d{k B T) 2 ' 1 ' 

The obtained results have been presented in figure [3] (B). On that basis, we have found the destructive ef- 
fects of the depairing electron correlations on the value of the specific heat jump at the critical temperature: 
[AC (T C )] M . =(U / [AC CZb)]„. =0 . 3 = 1.88. 

The high value of the electron-phonon coupling constant in the PtH compound means that the values of the 
determined thermodynamic functions will differ from the expectations of the BCS theory [33j . The most convenient 
way to estimate the differences between the results of the Eliashberg theory and the results obtained in the framework 
of the BCS model is to calculate the dimcnsionless ratios: 



T C C N (T C ) AC(r c ) 

Rh = i%(o) ' and Ro = Wpyj' ( } 

The BCS model predicts that the values of Rh and Rc are universal constants for all the superconductors and 
they are equal to 0.168 and 1.43, respectively [331 ]. In the case of the PtH compound, the following results have 
been obtained: R H g (0.145,0.156) and R c g (2.07,1.97), for n* g (0.1,0.3). On that basis, it is clear that the 
superconducting state in the PtH compound cannot be correctly described by the BCS model, regardless of the 
physical value of the parameter /x* . 



III. THE ELIASHBERG FORMALISM ON THE REAL AXIS 



A. The Eliashberg equations 

In the framework of the Eliashberg formalism the precise value of the energy gap at the Fermi level can be calculated 
on the basis of the course of the order parameter on the real axis (uj). In the cxamineded case, the form of the function 
A (uj) should be determined using the Eliashberg equation in the mixed representation (34|: 
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M 



<(uj + i8) = — ^2 [X(lj — iu> m ) - /j,*9 (lj c — \u>. 



and 

Z (lj + iS) 



m=-M 
f +00 



duj'a 2 F 







N ( LJ ) + / ( Lj' - LJ 



(11) 



(j) [uj — to + iS 



+00 



duj'a 2 F 



M 



l 71 V . 

H — > \(lj — il 

lj 8 ^— ' 

M rn=~M 



m 

LJ 



iV ( w ) + / f w + cj 
AT (V) +f(uj'-u 



\l 1 {to — lj') 2 Z 2 (lj — lj' + iS) — <j> 2 (cu — lj' + iS) 

LJ + LJ + id 



\J{lj + u'f Z 2 (lj + lj' + iS) - 4> 2 (lj + lj' + iS) 



(12) 



^LJ — LJ \ 



Z (lj — lj + id 



ITT 



"+00 ^ 

— / dui a 2 F ( lj 
^ Jo 



N (lj ) + f (lj +LJ 



\J ' (lj — lj' ) 2 Z 2 (lj — lj' + iS) — 4> 2 (lj — lj' + iS) 

(lj + U!^j Z ^LJ + LJ + iSj 
J '(lj + lj' 'f Z 2 (lj + lj' + iS) - <f> 2 (lj + lj' + iS) 



where the symbols N (lj) and / (lj) denote the functions of Bose-Einstein and Fermi-Dirac, respectively. 

The Eliashberg equations on the real axis have been solved using the numerical methods presented in the papers 
[IH, [1^. The convergent solutions have been obtained in the identical temperature range as for the Eliashberg 
equations on the imaginary axis. 



B. The energy gap 

Figure S] presents the order parameter form on the real axis. The frequency range, at which the Eliashberg function 
is defined has been chosen (lj £ (0, fi max )). We have considered the selected values of the temperature and the selected 
values of the Coulomb pseudopotential. 

It is easy to notice that the order parameter takes the complex values. The real part of the function A (lo) is 
used to determine the energy gap at the Fermi level, whereas the imaginary part of the order parameter determines 
the damping effects [37| ■ In the case of low frequencies, and thus the physical area of the energy gap, the following 
occurs: Re [A (lj)) Im[A(w)], which means small damping effects. For higher frequencies, the course of the real 
and imaginary parts of the order parameter is very complicated. However, for T = 2.5 K and //* = 0.1, the visible 
correlation between the shape of the order parameter and the Eliashberg function can be clearly seen. As it turns out, 
this correlation disappears with the increasing value of the temperature and the parameter (j,* . For the frequencies 
much higher than f2 ma x, the order parameter becomes the subject of the saturation. 

The energy gap at the Fermi level (2A) can be calculated using the equation: 

A(T)=Re[A(w = A(T))]. (13) 

From the physical point of view, the most important is the value of the energy gap at the lowest temperature (lo). 
In the considered case, the following result has been obtained: 2A (0) G (11.21,5.78) meV. Hence, the dimensionless 
ratio i?A = 2A (0) /ksTc changes in the range of the values from 4.25 to 3.98. It should be noted that the BCS 
theory provides the universal value of Ra for all the superconductors: [-/?a]bcs = 3.53 [H[. The result obtained for 
[i?A]p t fj confirms the fact that the properties of the superconducting state, can not be properly determined in the 
framework of the BCS model. 



IV. SUMMARY 



The study has characterized the thermodynamic properties of the superconducting state in the PtH compound 
under the pressure at 76 GPa. Due to the high value of the electron-phonon coupling constant, the calculations have 
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FIG. 4: The dependence of the order parameter on the frequency for the selected values of the temperature and the Coulomb 
pseudopotential. Additionally, the form of the rescaled Eliashberg function has been plotted. 



been carried in the framework of the Eliashberg formalism. A wide range of the Coulomb pseudopotential values has 
been chosen: /i* e (0.1,0.3). 

It has been found that the critical temperature is relatively high in the whole range of the Coulomb pseudopotential: 
Tc G (30.6, 16.8) K. The values of the remaining thermodynamic parameters differ significantly from the expectations 
of the BCS theory. In particular, the dimensionless ratios of the characteristic values of the thermodynamic functions 
arc respectively equal to: R H G (0.145,0.156), R c G (2.07, 1.97), and R A G (4.25,3.98). 
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